Fundamental scaling laws of on-off intermittency in a stochastically driven dissipative 

pattern forming system 
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Noise driven electroconvection in sandwich cells of nematic liquid crystals exhibits on-off inter- 
mittent behaviour at the onset of the instability. We study laser scattering of convection rolls to 
characterize the wavelengths and the trajectories of the stochastic amplitudes of the intermittent 
structures. The pattern wavelengths and the statistics of these trajectories are in quantitative agree- 



ment with simulations of the linearized electrohydrodynamic equations. The fundamental r 



-3/2 

distribution law for the durations r of laminar phases as well as the power law of the amplitude 
distribution of intermittent bursts are confirmed in the experiments. Power spectral densities of the 
experimental and numerically simulated trajectories are discussed. 

PACS numbers: 05.40.-a 61.30.-v 47.20.-k 47.54.+r 



I. INTRODUCTION 



A. Intermittency 



Intermittency is a prominent phenomenon observed in 
a large variety of nonlinear dynamical systems. The 'clas- 
sical' examples for intermittent behaviour, the so-called 
Pomeau-Manneville types I-III Jl|, || , can be found in de- 
terministic systems where upon a certain change of a con- 
trol parameter a fixed point of the system (corresponding 
to a periodic trajectory) becomes unstable. One of the 
characteristic features for the distinction of the different 
types of intermittency is the statistics of the duration 
of the quasi-periodic ('laminar') phases which are irreg- 
ularly interrupted by chaotic parts of the trajectory. 

A fundamentally different type of intermittent be- 
haviour has been observed in coupled chaotic oscillators 
j^, U ^). This phenomenon can be found in dynami- 
cal systems at the stability threshold when a stochastic 
or chaotic process couples multiplicatively to the system 
variables. The term on-off intermittency has been coined 
for this phenomenon. In systems which exhibit this type 
of intermittency, there is no sharp transition from an 
equilibrium quiescent state into an active state but in- 
termittent behaviour occurs for a range of values of the 
control parameter, and the system has to be character- 
ized by a statistical description. It resides in a ground 
(off) state during quiescent or laminar periods, which 
are interrupted by bursts of large scale excursions of the 
system variables into the on-state. As well as the other 
types, on-off intermittency is characterized by fundamen- 
tal statistical properties of the intermittent process which 
have been extensively studied during recentyears, exper- 
imentally as well as theoreticall y |rf , [^, || || |Io| , [TT] |l2] , 

H jg [10, 0, 0, ^ g |f gfp. a 

statistical analysis reveals characteristic asymptotic laws 
that describe the universal behaviour of such systems. 
It has been shown that the distribution of the durations 
r of the laminar phases in on-off intermittency follow a 
characteristic power law with exponent -3/2 Q in the 



vicinity of the instability threshold. Other fundamental 
scaling laws have been predicted for the distribution of 
the amplitudes of the bursts & 0j2(l_p7|and the power 
spectrum of the trajectories pi], |22j, |23|, |26[|. 

Experimental evidence for on-off intermittent behav- 
iour has been reported in a number of very different phys- 
ical systems. A simple experimental realization can be 
achieved in coupled oscillator circuits ||, other systems 
described in literature involve a gas discharge plasm a M 
and a ferromagnetic resonance spin wave experiment [p.0[ . 
While the fundamental validity of the asymptotic scaling 
laws is established theoretically, it is not easy to con- 
firm this prediction in the experiments. In the spin wave 
system | fTo| , a power law in the distribution of laminar 
phases has been reported over just slightly more than one 
order of magnitude in t. In the gas discharge experiment 
^| , the power law behaviour is covered by an exponential 
function. 

Among the experimental situations where on-off inter- 
mittent behaviour has been unambiguously detected is 
electrohydrodynamic convection (EHC) in nematic liquid 
crystals driven by multiplicative noise jn|. This system 
turns out to be particularly well suited for an experimen- 
tal characterization. It represents the first spatially ex- 
tended dissipative system where on-off intermittency has 
been identified. Compared to other reported systems, the 
EHC experiment shows an additional spatial periodicity 
of the intermittent bursts where a wavelength selection 
process is involved. The access to the control param- 
eters and observation of trajectories is straightforward. 
The physical mechanisms are well understood. Many in- 
volved material parameters are accessible in independent 
experiments. The validity of the asymptotic scaling law 
for the duration r of laminar phases has been confirmed 
experimentally fllj] . 

In this study, a modified optical setup is used in order 
to record pattern wavelengths and amplitudes with high 
sampling rates. The trajectories of pattern amplitudes 
are extracted from the laser scattering profile produced 
by the nematic cell. In addition, a simulation of the lin- 
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earized dynamic equations is presented. Trajectories ob- 
tained in the simulations are compared to the experimen- 
tal results to test the validity of a linearized treatment of 
the system dynamics [E8| near the instability threshold. 



B. Nematic Electroconvection 

Nematic EHC @ fo[ ||, |§ §| |§ represents a stan- 
dard system of dissipative pattern formation, its funda- 
mental features are well understood today p5| . The in- 
stability is driven by interactions of an external electric 
field with electric charges, present as either impurities or 
dopants in the nematic material. 

The experimental geometry is sketched in Fig. ^. The 
ground state with the nematic director (optic axis) uni- 
form along an easy axis in the cell plane is achieved by 
proper surface treatment of the cell plates. An electric 
field is applied normal to the cell plane. Nematic mate- 
rial with negative dielectric anisotropy Ae is chosen to 
prevent the splay Freedericksz instability. In the electric 
field, the dielectric torque stabilizes the ground state, but 
small thermal fluctuations of the director field in connec- 
tion with an anisotropic conductivity of the nematic gen- 
erate a periodic space charge modulation in the cell plane. 
The interaction of these space charges with the electric 
field lead to convective flow which in turn generates a 
destabilizing viscous torque on the director field. Above 
a threshold voltage U c , this torque exceeds the stabiliz- 
ing dielectric and elastic torques and a periodic pattern 
of convection rolls and corresponding director deflections 
is formed. In the simplest case, the wave vector of the 
first instability is along the director (normal rolls). The 
roll structure appears in optical transmission as an array 
of parallel stripes in the cell plane (Fig. []]). 

Electroconvection is conventionally driven with a pe- 
riodic AC voltage (to avoid sample degradation in DC 
fields). For the understanding of the mechanism of pat- 
tern formation it is essential to note that director and 
charge fields respond on different time scales to the al- 
ternation of the electric field. The symmetry of the dy- 
namic equations requires that their time dependence with 
respect to a periodic driving field is antisymmetric. This 
is reflected in the existence of two different types of pat- 
terns. In the low frequency 'conduction' regime, charge 
relaxation is fast compared to the AC frequency. The 
charge density alters its sign with the applied field, the 
sign of the director deflection is preserved. The voltage 
for the onset of conduction rolls increases monotonously 
with frequency. At the cut-off frequency / c , the threshold 
voltage curve intersects that of the 'dielectric' patterns. 
In the dielectric regime above / c , the director deflections 
alternate with the field while the charge density modula- 
tion retains its sign. Figure || shows the stability diagram 
for the system studied in this paper. Figure y(a) gives 
the onset voltage for the first instability, towards nor- 
mal rolls. The corresponding wave numbers are shown in 
Figure |(b). 



When a stochastic excitation scheme is used where the 
driving field has no deterministic component (as e.g. the 
dichotomous Markov process that is considered in this 
paper) the system does not exhibit a sharp transition 
from the quiescent to the convective state upon varia- 
tion of the control parameter but shows intermittent be- 
haviour. Two different regimes are found which have 
many features in common with the corresponding con- 
ductive and dielectric regimes in the deterministic case. 
At onset of the instability, one of the system's character- 
istic times (director or charge relaxation) becomes com- 
parable to the characteristic time of the noise T st0 ch ■ A 
considerable qualitative change of onset and appearance 
of the convection patterns J 36|, [57| is observed. A typical 
snapshot is shown in Fig. |3. 

Nevertheless, there was little interest until recently in 
the quantitative statistical interpretation of the struc- 
tures at the instability threshold. Main focus of research 
has been directed to the study of superimposed deter- 
ministic and stochastic driving fields and the construc- 
tion of pattern state diagrams. The empirical concept of 
a threshold voltage has been applied in previous exper- 
imental investigations of noise driven EHC |36|, |3^, |3^, 
51 ^2) and statistical methods have been applied to 
test various stability criteria ||, ||J, ||, Q. The 
largest Lyapunov exponent A of the trajectories of the 
system variables, which is analytically known p8|, pro- 
vides a quantitative measure for the instability threshold, 
but the crucial problem is the experimental determina- 
tion of the threshold in a system with limited dynamical 
range and additive noise. The statistical analysis of the 
intermittent trajectories provides an experimental tool to 
characterize the stability threshold, A = 

The observation and quantitative characterization of 
the dissipative patterns under stochastic excitation is 
achieved in this study by exploring the phase grating 
formed by the spatially periodic deformed director field. 
The laser scattering profile resulting from the nematic 
convection patterns is recorded. A two-dimensional (2D) 
detector allows to observe the complete scattering profile 
of the transmitted laser light and to study the wave vec- 
tor orientation and wavelength selection process with a 
sampling rate of 25 s , and alternatively a photo diode 
positioned at the scattering reflex enables us to record 
the trajectory of the dominant mode with faster sam- 
pling rate and higher intensity resolution. 



II. SAMPLE PREPARATION AND 
EXPERIMENTS 

A. Sample preparation 

We use the nematic mixture Mischung 5 (Halle) 
which consists of four disubstituted alkyl(oxy)phcnyl- 
alkyloxybenzoates. This material has been used in pre- 
vious EHC experiments JO], |37], [l7|, |||. The mixture 
is nematic at room temperature, its clearing point is 
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70.5 °C. The commercial sandwich cells (LINK AM) used 
in the experiments have cell gaps of 25.8 /im. The glass 
plates are covered with transparent Indium-Tin-Oxide 
electrodes (5 mmx5 mm), they are polyimide coated and 
antiparallelly rubbed for planar surface alignment. The 
temperature of the samples is controlled by a Linkam 
heating stage with an accuracy of 0.1 K. The sample 
temperature has been set to 32°C in all experiments. 



B. Material Parameters 

All relevant material parameters for the simulation of 
the electrohydrodynamic equations except for some vis- 
cosities have been measured in independent experiments 
p9| , |50| . The conductivity of the nematic samples dif- 
fers by about 20% between individual cells, values for a 
given cell are almost constant in time. In order to pre- 
vent long-term trends of the conductivity, we reheat the 
material into the liquid phase between subsequent runs of 
the experiment, similar procedures have been proposed 
in |5lJ ^2[ . All measurements presented in the diagrams 
of this paper have been performed consecutively with the 
same cell, in order to obtain quantitatively comparable 
results for the different statistical investigations. Differ- 
ences in conductivity between individual cells may lead 
to variations of the respective thresholds but do not in- 
fluence the statistical characteristics. 

In order to complete the parameter set for the nu- 
merical simulations, we fitted the threshold curves and 
critical wave numbers for deterministic square wave exci- 
tation where the experiment yields sharp thresholds to- 
wards the first instability. With a fixed parameters set 
(Table I) , good agreement between numerical results and 
experimental data is achieved across the whole frequency 
range investigated (Fig. ||). The same parameter set is 
used thereafter for the calculation of the Lyapunov ex- 
ponents and stochastic trajectories. 

C. Excitation 

The waveform of the driving electrical voltage is syn- 
thesized by a computer and subsequently amplified. We 
generate the driving voltage curve with a sampling rate 
of 10 kHz. Although the spectrum of this waveform is in 
principle discrete, we can consider it as quasi-continuous 
in the frequency range relevant in the experiments (below 
500 Hz). Within this study, we have used a few special 
waveforms, which are detailed in the following, Figure ^ 
visualizes their signal shapes and shows the respective 
power spectra. 

Periodic square wave excitation is used to construct 
the stability diagram for deterministic driving. We ex- 
ploit the sharp threshold towards the first instability for 
the adjustment of system parameters and test of the long 
term stability of the samples. For example, the cut- 
off frequency is a sensitive measure to reveal even small 



changes of the sample conductivity. 

The Dichotomous Markov process (DMP) is the sto- 
chastic waveform used in all experiments presented be- 
low. It is characterized by random jumps of the elec- 
tric field between the two values +E and — E, with an 
average jump rate a. The time intervals between con- 
secutively jumps are distributed as Ati = — 1/a lna^, 
where xi S (0, 1] is a uniformly distributed random num- 
ber. In analogy to a deterministic square wave excita- 
tion with frequency v, we define the 'average frequency' 
v = a/2. The DMP power spectrum is Lorentzian with 
its maximum at frequenzy zero and a half width related 
to the jump rate by a /it (Fig. ||(b)). An important fea- 
ture that facilitates the numerical calculations is that all 
terms quadratic in the electric field are time independent 
in the DMP. An analytical calculation of the sample sta- 
bility threshold and the critical wave numbers have been 
performed [^8| . Technically, identical realizations of the 
stochastic process can be reproduced with the computer. 
This allows us to use identical noise sequences for exper- 
iment and simulations when details of the trajectories of 
pattern amplitudes are of interest. 

Other stochastic processes have been tested in addi- 
tional experiments. While the DMP randomizes the 
phases of a periodic square wave, another stochastic pro- 
cess can be synthesized that randomizes the amplitudes 
of the square wave but keeps the jumps equidistant. This 
randomization of the amplitude can be combined with 
a random phase of the jumps. Such processes do not 
complicate the numerical simulations, the electric field is 
piecewise constant and the equations of motion can be 
integrated straightforwardly. It turns out that the statis- 
tical properties of the trajectories and the derived scaling 
laws for these stochastic processes are in full accordance 
with those for the DMP. Therefore, we will discuss only 
the DMP results being representative for other types of 
stochastic excitation. 



D. Data acquisition 

The sample cell is irradiated by a He-Ne laser with 
wavelength Ai aser = 632.8 nm at normal incidence. The 
beam diameter is about 1 mm. The scattered light is 
monitored on a screen in a distance of maximum 1.3 m 
from the cell (Fig. |J). When the sample is in the ground 
state (zero electric field), only a weak background scat- 
tering Ao(t) is observed around the primary beam. Since 
we are interested in scattering from spatially deforma- 
tions of the director field, we correct the raw data with 
a time averaged intensity Aq by A(t) = \A Taw (t) — Aq\. 
This correction is marginal since the background signal 
is in general three orders of magnitude smaller than the 
amplitude of the scattering reflex at the position of the 
photo diode. 

From the two-dimensional scattering images, wave- 
lengths and orientations of the patterns can in principle 
be continuously extracted. However, because of band- 
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width limitations in signal processing we use two different 
equipments to record the scattering data. 

A commercial video camera is employed to take 2D 
scattering images with an acquisition rate of 25 frames 
per second and 8 bit intensity resolution. This enables 
us to study the evolution of the mode spectrum and to 
access the full wave vector information, although time 
and intensity resolution are limited. 

The 2D images show that at given frequency, the pat- 
tern is dominated by a single mode (see Fig. ^) with fixed 
wave vector but varying amplitude. Any low order scat- 
tering reflex of that mode is representative for the mo- 
mentary pattern, and it is sufficient to record only the 
scattering intensity at a fixed position (here, we use one 
of the two symmetric second order reflexes). For that 
purpose we employ a photo diode with aperture 7 mm 2 
adjusted to the reflex of interest. 

The photo diode signal is digitized and the trajectory 
is stored in a computer. We use a fast (0.2 ms time 
resolution) 12 bit AD concerter for measurements with 
high time resolution as for example the determination 
of laminar phases, and alternatively a slow (0.1 s) 24 
bit ADC for accurate amplitude measurements over large 
dynamic ranges. 



III. OPTICS 

The electroconvection rolls in the nematic material 
produce spatial periodic director deflections 



without deflection, creating a phase grating 



<p(x, z, t) = tp(t) sin(fe x x) sin(fc 2 z) 



(1) 



which lead to a modulation of the effective refractive in- 
dex n c g for the transmitted extraordinary beam, 



n e s(P(x, z)) 



\J n 2 cos 2 P(x, z) + n 2 sin 2 (3{x 1 z) 



, (2) 



where n Q ,n e are material parameters and (3 is associ- 
ated with the angle between the light ray and an opti- 
cal axis []53|]. A phase and amplitude grating is formed 
in the cell. There has been some discrepancy in litera- 
ture about the usage of refractive or ray index in these 
calculations. For the small fluctuations considered here, 
both approximations lead to similar results. In order 
to establish a relation between the experimental obser- 
vations and the results of the simulations of the direc- 
tor and charge fields, it is necessary to calculate the 
diffraction efficiency of a given periodic director field. We 
use the eigenfunctions tp{x,z,t) from the linear stability 
analysis. Using Format's principle it is possible to de- 
termine the optical path, the resulting phase difference 
4>(x, t) and the intensity profile of light passing the cell 
[§§ ||, ||, H, ||, || |6Cj. At small deflection angles, 
one can assume that individual light rays pass the cell 



laser 



(j>{x,t) 



1 
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"off dz 



laser JQ 
,2 
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(3) 



2n 2 



-tp(t) 2 cos 2 (k x x) sm 2 (k z z) 



Analytical integration over z along a straight path 
through the cell pl[ | yields a periodic phase modulation 
A<f)(x, t) with twice the wave number k x and a quadratic 
dependence on the director deflection 



Ac/>(x,t) = 



n e 7T 2 {n\ - n 2 ) 



ip(t) 2 cos 2 (k x x) 



2?1 2 fc z Al aS er 

oc ip(t) 2 cos 2 (k x x) oc <p(t) 2 cos(2 k x x). (4) 

In this approximation, the intensity modulation with 
k x due to focussing effects of the inhomogeneous refrac- 
tive index profile is neglected and only even order reflexes 
appear in the scattering image. This is in agreement with 
the experimental observations. The second order reflex 
dominates the small amplitude patterns, and with in- 
creasing amplitude of the director deflections tp(t), higher 
order reflexes can be observed. We note that the ampli- 
tude grating (which produces the well known shadow- 
graph images in conventional orthoscopic microscopy) is 
effective as well, it is most prominently reflected in the 
weak first order reflexes (see Fig. |^). 

The relation k x — 2nsmO m /(mXi ascl -) connects the 
scattering angle m of the mth order reflex with the wave 
number k x (8 m ). The scattered light intensity A(t) of the 
phase grating at the second order reflex is related to the 
square of the Bessel function of first kind, Ji((f>max(t)) , 
with the amplitude of the phase grating ma x in the ar- 
gument Hi, II H H, H ||. For small director de- 
flections tp{t), the intensity at the second order reflex is 
proportional to the fourth power of tp 



A(t) oc J 2 1 



(*)) oc <t>l ax (t) « v\t). 



(5) 



By numerical integration of the nonlinear Euler- 
Lagrange equations we have calculated the actual propa- 
gation of light beams through the sample. The numerical 
integration allows for the deflection of light and thus for 
intensity modulations by focussing effects. These simu- 
lations confirm the relation A(t) oc tp(t) A even for large 
director deflections tp(t) < 1. 

The numerical simulation of EHC patterns yields ip = 
d x (p, (see below) which represents the director deflection 
for a given mode. For comparison of simulated tp and 
experimental intensities we use the relation 



A th (t) = constat). 



(6) 



We emphasize that we compare the experimental and 
simulated pattern amplitudes only on a relative level, 
since no efforts have been made to relate calculated ab- 
solute scattering intensities to voltages measured by the 
photosensor. Since the simulations use a linear model, 
absolute scaling of pattern amplitudes is not relevant for 
the fundamental statistical properties of the system. 
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IV. NUMERICAL SIMULATIONS 

Basic Equations 

The basic equations for the charge and director fields 
and the method for analytical derivation of the Lyapunov 
exponents in the electrohydrodynamic instability have 
been described in detail in pq] . 

The important quantities describing the system dy- 
namics are the charge density q and the amplitude of the 
director deflection tp (see Fig. Q). A standard technique 
to describe the time evolution of small amplitude peri- 
odic patterns is the use of linearized differential equations 
and a two-dimensional mode ansatz for the two involved 
quantities, 

ip(x,z,t) = ip(t) sm(k x x) sin(fc z z), (7) 
q(x,z,t) — q(t) cos(k x x) sin(fc z z). (8) 

In the relevant parameter range, the pattern is spa- 
tially periodic along the director easy axis (x direction) 
This result of the linear stability analysis is in agree- 
ment with the experimental observations. Therefore, 
we consider only the stability of modes with the wave 
vector parallel to the x axis. The (stress free) bound- 
ary conditions for the director field at the glass plates 
(p(z — 0) = (p[z = d) = are satisfied by k z = n/d. 
For convenience, the director deflection (p is substituted 
by the curvature ip = d x (p. A system of two ordinary 
differential equations is derived from the torque balance, 
Navier-Stokes and Maxwell equations. After lineariza- 
tion we obtain an ordinary differential equation system 
in time for the vector z (t) — (q (t) , ip (t)) T , 

-z(t)--( 1/T « a * E{t) Az(t) (9) 
dt [ ) ~ \aE(t) A 1 -A 2 E(t) 2 J zW ' W 

where Ai, A2, T q , <jh and a are parameters related to the 
viscous, elastic and electric properties of the liquid crys- 
tal as well as to the wavelength k x of the modes ||^] . The 
electric field amplitude E(t) =U(t)/d corresponds to the 
excitation voltage U(i) . In case of piecewise constant ex- 
citation (as the deterministic square wave and DMP de- 
scribed above), E{t) assume only two values ±E, and all 
elements of the matrix are constant between consecutive 
jumps. Within a time interval Ati of constant excita- 
tion field, integration of the differential equation gives a 
solution in the form of a sum of two exponentials. The 
solution Zj from the ith integration step is taken as the 
starting value for the (i + l)st step. The complete tra- 
jectory for given k x , E and set of (random) times Ati 
for jumps of the excitation voltage is calculated with an 
initial vector z$ = {q,ip)\t=o- At the discrete jump times 
t n , the solution is 



z(t n ) = T<»>(At n )-..T<«*>(Ati)z(0), (10) 
AU = U — U-i ; Si = 8ignE(ti > t > U-i), 

n 

r n = J] T^(Ati), (11) 

»=1 

where T^ Si '(Ati) is the 2x2 time evolution matrix for the 
ith interval. 

For a detailed statistical analysis oiip(t), in particular 
for the calculation of power spectral densities (PSD), the 
trajectory can be filled in the intervals between the jumps 
using the exact exponential solutions T^ Si \t — U-\). 

In the particular system studied here, the eigenval- 
ues Eigenval7^, (_E, k x ) are real and positive. For pe- 
riodic square wave driving, all intervals Ati are equal 
and the product is T n = (T+T~)™/ 2 . For the calcula- 
tion of the Lyapunov exponents it is sufficient to consider 
Eigenval (T + T _ ). This reproduces the well known re- 
sults of classical theory using Floquet methods |3^, |3j| . 

In the case of stochastic excitation all the Ati are dif- 
ferent and the calculation of the Lyapunov exponents 
leads to an infinite product of 2 x 2 random matrices |28| . 
This system yields two real Lyapunov exponents Ai > A2 
which are related to the eigenvalues of the product of 
stochastic matrices T n . In particular, the largest Lya- 
punov exponent (in the following denoted by A) is found 
from 



A= lim — ln{max(Eigenval7^)}. (12) 

rwoo t n 

For DMP excitation, the Lyapunov exponents can be 
obtained analytically |2£| . Figure [7| shows the Lyapunov 
exponent for the critical wave length calculated with the 
parameters specified above. In the linear model, |z| is 
growing to infinity (A > 0) or shrinking to zero (A < 0), 
depending on the value of the largest Lyapunov expo- 
nent. The selected wavelengths and threshold voltages 
for a given frequency and set of material parameters are 
calculated from the neutral curve. The wave number is 
varied and the minimum E c of the neutral curve provides 
the critical wave number k c . 

Because of the symmetries of Eq. (^), one of the sys- 
tem variables (q, ip) keeps its sign while the other variable 
has to change its sign with the polarity of the applied 
field. At periodic excitation, the system is synchronized 
with the applied field and the conduction (q alternating 
periodically) and dielectric (ip alternating periodically) 
regimes are distinguished. In case of DMP excitation, 
the support of one of the variables still preserves the sign 
and the two regimes can still be distinguished |^8| . In the 
following, we will discuss only the low frequency conduc- 
tion regime, Fig. || shows details of a simulated trajec- 
tory. The slow variable is the director deflection which is 
directly related to the measured quantity, the scattering 
intensity. Due to the coupling of (q, ip), both the director 
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deflection and the envelope of the charge density curves V. RESULTS AND DISCUSSION 

show the same long term (t > 1/v) characteristics. 

A. Pattern images 



Boundaries 

In the linear description, the trajectories tend to infin- 
ity or zero for all values of A ^ 0. However, a realistic 
description of the experiments has to consider boundaries 
for z. Nonlinearitics in the system limit the excursion of 
the system variables z to large values while additive noise 
prevents their unlimited decay. In order to compare the 
results of simulations and experiments, we introduce lim- 
its for the director deflection (curvature tp) by clipping 
each time step 

max • 

Ipmin if 1p(U) < -0min, (13) 

ip(ti) else. 

Because of the linearity of the ODE system (Eq. (Bh), 
only the ratio of the upper and lower thresholds is impor- 
tant. A constant factor in the amplitudes is irrelevant for 
the statistical properties and scaling laws. Here, we as- 
sume that the dynamic range is two orders of magnitude, 
5 x 10 -3 < ip(U) < 0.5. This dynamical range reflects ap- 
proximately the situation of a thermal background stim- 
ulation {ip 2 ) 1/2 ~ ii/K n kl ~ 5 mrad @ |f, || and 
an upper limit of 0.5 rad. For negative Lyapunov ex- 
ponents, the lower boundary is necessary to prevent the 
unlimited decay of ip(t). The choice of the value of the 
lower boundary has strong effect to the number of burst 
per unit time (see Fig. ||), but only small effects on the 
fundamental statistical laws (see below). 

The assumption of a well defined lower limit is of course 
artificial and cannot describe the actual experimental be- 
haviour for very small pattern amplitudes. A more realis- 
tic assumption considers low amplitude additive noise in 
the variable ip(ti). We have studied this case by adding 
Gaussian random numbers with a given variance D Q 
at each time step. For zero electric field, Eqs. (||) de- 
couple and -0 describes an Ornstein-Uhlenbeck process 
(OUP) which is characterized only by D and the expo- 
nentially decaying autocorrelation with the correlation 
time roup oc 1/A related to the Lyapunov exponent at 
zero voltage (see Fig. Q). One consequence of such addi- 
tive noise is that the simulated trajectories for the slow 
variable ip(t) can change sign (see Fig. [To]), in contrast to 
Eq. (|l^) . This is not relevant for the comparison with the 
scattering experiment which is insensitive for the spatial 
phase of the pattern. If the variance D is chosen such 
that the mean square amplitude of tp in absence of elec- 
tric fields gives (ip 2 ) 1 / 2 = -0min, the statistical properties 
of | "01 are qualitatively identical with the simulations us- 
ing conditions (|l3|). 



The conventional technique to observe pattern forma- 
tion in EHC is the shadowgraph method Jm] in combina- 
tion with a transmission microscope. An instant picture 
of a pattern burst, recorded by means of orthoscopic mi- 
croscopy, is shown in Fig. ||(a). The dynamics of this 
pattern can be visualized best when the intensity in a 
cross section perpendicular to the rolls is scanned and 
a spatio-temporal plot as in Fig. |3|(b) is analyzed. Some 
bursts, in particular those which appear in fast sequence, 
are correlated in their spatial phase. After long laminar 
phases, however, there is in general no spatial correlation 
remaining between the bursts. For stationary rolls, this 
loss of correlation is a consequence of additive noise in 
the system, it triggers new modes whenever the pattern 
amplitude reaches thermal noise level. In such cases, the 
new fluctuation mode has an arbitrary spatial phase with 
respect to the previous convection pattern. On the other 
hand, there are sequences of bursts in Fig. ||(b) which 
preserve spatial correlation. They can be attributed to 
amplifications of the same mode which has disappeared 
in the optical image but has not reached the level of ad- 
ditive noise during the intermediate laminar phase. 

B. Scattering images 

Figure ^(a) shows a snapshot of the scattering image 
originating from a burst in a DMP driven d=25.8 /xm 
thick nematic cell, taken with the 2D camera detector. 
The image shows the primary beam at 9 = and dif- 
fuse scattering spots from the periodic spatial director 
modulation in the cell. The scattering reflexes concen- 
trate on the k x axis, i.e. normal electroconvection rolls 
are observed. A wave number k x of the pattern pro- 
ducing this image of 0.2 /jm" 1 is derived from the spot 
positions, the director period is Xdir = 2Tr/k x = 32 /im. 
Below the scattering image, in Fig|^(b), the time depen- 
dence of the intensity profile taken along the horizontal 
symmetry axis of the profile (k y = 0) is plotted. For bet- 
ter reproduction, the image is plotted with inverse gray 
scale, dark spots correspond to high amplitudes of scat- 
tered light and consequently to large director modula- 
tions (bursts). The bursts are characterized by a narrow 
wavelength band and the reflexes remain approximately 
at the same positions, i.e. in all bursts the patterns have 
nearly the same wavelengths. The wavelength selection 
process itself is not observable in the images, it is obvi- 
ously fast compared to the video rate of 25 Hz and occurs 
below the level of optical sensitivity of our camera. The 
information about the spatial phase of the patterns gets 
lost in the scattering image, therefore we cannot deter- 
mine from the scattering images whether the convection 
rolls of subsequent bursts appear spatially phase corre- 
lated or not. The essential information taken from the 2D 
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images is that the wavelength of the noise driven patterns 
is constant and the scattering image consists of reflexes 
at fixed scattering angles with varying amplitudes. Thus, 
we achieve a considerable data reduction by restricting 
to the strongest scattering reflex. In particular, the de- 
tector is placed at the second order scattering reflex of 
the most critical wave number, indicated by the arrow in 
Fig. |(b). 



C. Trajectories 



The intensity of scattered light at the position of the 
second order scattering reflex of the first unstable mode 
is shown in Fig. O. The curves have been digitized by 
means of the 24 bit AD converter. The random elec- 
tric excitation field E t uses identical realizations of the 
stochastic process for all three trajectories, with different 
amplitudes E. Concerning the effects of the multiplica- 
tive noise in the system, details of the three curves can 
be directly compared. The characteristic frequency v was 
80 Hz, corresponding to an average of 160 jumps/s of the 



sign of E t . Fig. 



ll|(a) shows the raw signal from the 

The con- 



detector for a voltage below the critical U° xp 
tributions of (additive) background fluctuations to the 
detector signal are of the order of 10 -3 V around a con- 
stant offset A = 3 x 10~ 3 V. Bursts of the spatially 
periodic pattern that exceed the background noise level 
occur infrequently. From arguments discussed in the next 
sections we conclude that the excitation voltage is below 
the stability threshold (defined by the Lyapunov expo- 
nent A = 0). In Fig. [Tl](b), the voltage is approximately 
equal to the critical voltage (A = 0). The characteristic 
feature of the trajectory is its mirror symmetry of high 
and low amplitude excursions of the scattering intensity 
A(t) on logarithmic scale. Rise and decay sections of the 
graph are symmetric. We note that in a linear presen- 
tation of the same plot A(t), the high amplitudes will 
appear as prominent bursts out of the background level, 
and the typical intermittent behavior is acknowledged. 
At higher voltages (Fig. [ll](c)), the trajectory will be 
predominantly in the saturation region ('on' state), inter- 
rupted by short laminar phases. These intrinsic symme- 
tries reflect theoretical predictions for on-off intermittent 
behaviour p9| . 

Amplitudes in Figs. [n](b,c) have been corrected for 
the background intensity, A = \A Iam — Aq\. This cor- 
rection attempts to separate the constant (stray light) 
background signal and additive (thermal) noise in the 
trajectories. Since these contributions are comparably 
small, the correction effects only the low amplitude sec- 
tions of the trajectories. It faciliates the identification of 
laminar phases and enables us to visualize the symmetry 
of burst and quiescent periods in the logarithmic scale. 



D. Distribution of Laminar Phases 

The statistics of experimental and simulated trajec- 
tories can be compared on a quantitative level. In the 
experiment, the distribution of laminar phases is calcu- 
lated by introducing an arbitrary threshold intensity A on , 
and the durations of periods where the intensity curve 
stays completely below that threshold are determined. 
A threshold A on = 0.05 V has been chosen here, it cor- 
responds to the geometrical average of lower and upper 
bound of the photo sensor signal. As expected in on-off 
intermittency, the choice of the actual threshold value is 
not critical. Fig. [l2] (left) depicts the distribution of the 
laminar phase durations extracted from trajectories for 
DMP excitation with three different voltages, the same 
values as in Fig. [n]. The time axis is scaled with the jump 
rate a of the DMP. These distributions were extracted 
from 4200 sec runs of the experiment, each trajectory has 
been recorded with a sampling rate of 1000 s _1 (with the 
fast AD converter), so that a range of 6 orders of magni- 
tude in r is resolved. The dash-dot line indicates a t~ 3 / 2 
dependence, which is predicted theoretically exactly at 
the sample stability threshold, A = 0, when no additive 
noise is present. In the short time range, for 2vt < I, 
the curve deviates from this predicted fundamental de- 
pendence, because one approaches the time scale of the 
driving DMP process and specific details of the driving 
process become important. In the long time limit of the 
curve, the power law behaviour of the experimental tra- 
jectories breaks down mainly because of the lower bound- 
ary (additive noise level) for the system variables. For 
long periods at least one of the variables (q, ip) reaches 
a (thermal) noise level which prevents excursions of the 
system variables to values much below A on , trajectories 
are practically reflected there. This lead to faster injec- 
tions of the next burst above A on , and thus to a lower 
probability of long duration periods. The flat shoulder 
indicated in Fig. [l^ is the outcome of this effect Q . 

Figure [T^ (right) shows the results of the corresponding 
simulations. Limits to the system parameters were set as 
given in Eq. (|l|) . The voltages used in the simulation are 
U' h =I4.2 V (corresponding to the critical value A = 0), 
and U' h ±0.7 V chosen in the vicinity of this threshold. 
The critical voltage in the simulated curves is derived 
from the analytically calculated Lyapunov exponent [ p8[ . 
It is in perfect agreement with the numerical simulation 
of the trajectories in absence of upper and lower limits. 
In the experiment, it has been proposed earlier that a 
reasonable definition of the critical voltage can be found 
when the statistical distribution of the laminar phases 
duration is analyzed. Thus we assume that the experi- 
mental critical driving voltage U® xp is reached when the 
distribution of laminar phases durations is most adapted 
to a t -3 / 2 dependence Knowledge of the relation 

(||) between director deflection amplitudes and scatter- 
ing intensities is not necessary for the determination of 
the distribution of laminar phases. On the other hand, 
the laminar phase distribution is not the most sensitive 
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criterion for the determination of the sample stability 
threshold as will be shown in the next section. 



E. Distribution of Pattern Amplitudes 

Another fundamental prediction for the statistics of 
on-off intermittent processes is the distribution of ampli- 
tudes in the trajectory. It has been shown ||, |I| ^6], 27j 
that the distribution of the amplitudes of the intermit- 
tent variable A should follow a power law 



p(A) oc A 



-l+JJ 



(14) 



in the vicinity of the threshold A = 0. The parameter 
rj oc A vanishes at the threshold of the instability. 

A statistical analysis of the recorded trajectories re- 
quires knowledge of the quantitative relation (|^) between 
the measured scattering intensity A and the amplitude of 
the pattern ip. One can easily show that a similar power 
law as for the amplitude distributions of the intermittent 
variable ip holds also for quantities A{p) that depend on 
the mth power of ip. When p(p) oc (p~ 1+r> and A oc p m 
then 



p(A) oc A- 1+ri ' m . 



(15) 



The hyperbolic dependence at the threshold A = is 
preserved. Our observable A(t), the scattered intensity 
in the 2nd order reflex, is related to the director deflec- 
tion amplitude in first approximation by Eq. (^) , and this 
relation provides another opportunity to determine the 
critical voltage U c . Fig. [l3| shows the distribution of scat- 
tering amplitudes A(t) for zero driving voltage and four 
representative voltages in the vicinity of the stochastic 
threshold U c . A power law can be fitted to the middle 
part of all these distributions. For low scattering ampli- 
tudes, the curve deviates from the fit because of the su- 
perimposed background scattering. For large amplitudes, 
the power law breaks down because of the saturation of 
the system variables and because the optical character- 
istics, Eq. (||), is not valid for large director deflections. 
In the numerically simulated trajectories, Fig. [l4[ where 
hard boundaries, cf. Eq. ( ^3[ ) have been used, these ef- 
fects are condensed in the edges of p(A). A th is computed 
from <p by means of Eq. (^) . 

In Fig. ||(b), the amplitude for U=12.9 V is closest 
to the exponent -1 in the power law, therefore we as- 
sign this voltage to the experimental threshold voltage 
U° xp . This value is in good coincidence with the criti- 
cal voltage found from the best fit of the laminar phase 
distributions to a r~ 3 / 2 law. Both statistical definitions 
of the experimental threshold voltage agree consistently, 
and the numerical simulations confirm the equivalence of 
the thresholds determined from the distributions of am- 
plitudes and laminar phase durations. 

At voltages near the threshold, the power law still 
holds and the exponents — 1 + rj can extracted. Typi- 
cal amplitude distributions are depicted in Figs. 0(c,d). 



The exponents extracted from experimental data as well 
as from simulated trajectories are shown in Fig. In 
accordance with the natural scales of the system, the 
axis was normalized with the critical voltage to obtain a 
control parameter e. The optical relation (^) has been 
applied to retrieve pattern amplitudes A which can be 
compared to the simulation from experimental scatter- 
ing intensities A. The good agreement in the slopes of 
the 77(e) curves at e > justifies the application of the 
optical relation Eq. (H). 



F. Power Spectral Density 

Finally, we discuss the power spectrum of the intermit- 
tent process. The main general theoretical prediction for 
the power spectral density (PSD) for A = is a square 
root frequency dependence ||l], [2^, ^6) in a certain 
frequency range. If any quantity with a power law de- 
pendence on the intermittent variable is observed, this 
prediction is equally valid ]2^ |. The PSD predicted for 
a process driven with multiplicative noise is qualitatively 
different from a process where noise couples additively to 
the system variables. 

In the range of very small frequencies, a constant PSD 
is expected because any time correlation in the system 
variables is destroyed by additive noise and the limited 
dynamical range of (q,ip). In the high frequency tail, a 
PSD oc u)~ 2 dependence is expected, similar to that of 
a simple stochastic process with exponentially decaying 
autocorrelation function. The relation for the high fre- 
quency limit u> — > 00 can be obtained analytically (from 
Eq. (50) in H) 



PSD oc ^ 



const if u! < oj\, 
u;- 1 / 2 if uji < uj < u>2, 
uj 2 



(16) 



if UJ > U>2. 



The crossover frequencies ui\ , ui2 where the asymptotic 
exponent changes depend upon the Lyapunov exponent 
and specific properties of the additive and multiplicative 
noise. 

Fig. [l^ displays the PSD of experimentally recorded 
trajectories A(t) for three voltages (same as in Fig. fl2| ). 
It was obtained from Fourier transformation of 4.2 xlO 6 
data points of the optical trajectories in an interval of 
4200 s. In the logarithmic plot, the curves have been 
smoothed by averaging the spectral energy over intervals 
Auj proportional to uj. 

At the threshold voltage U° xp (b), the uj' 1/2 and 
regions are well separated, and the existence of a constant 
PSD in the low frequency wing seems to be indicated. In 
the high frequency wing, influences of the mean frequency 
of the driving process, v = 80 Hz, are observable. 

The PSD obtained from the numerical simulation (see 
Fig. p"7| ) are not in agreement with the experiment, in 
particula the square root dependence is not reproduced. 
Only when the dynamical range is chosen unrealistically 
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large, we obtain a PSD oc uj^ 1 / 2 . The choice of a realis- 
tic lower bound (Eq. ([l3|)), i.e. additive noise with rea- 
sonable amplitude, destroys any long time correlations. 
Therefore, a simulation of the PSD with similar param- 
eters as in Figs. [l4| yields a pronounced constant low- 
frequency region. 

VI. SUMMARY 

On-off intermittency in stochastically driven electro- 
convection of nematic liquid crystals in the conductive 
regime has been investigated experimentally and by nu- 
merical simulations. Results have been presented for 
excitation with the dichotomous Markov process, but 
the resulting fundamental statistical behaviour is qual- 
itatively similar for many other types of stochastic exci- 
tations. 

Laser scattering has been used to determine the wave- 
lengths and time resolved pattern amplitudes. It has 
been shown experimentally, and confirmed in the simu- 
lation of the electro-hydrodynamic equations that under 
stochastic excitation the pattern selects its wavelength 
within a narrow band, therefore the intensity of scat- 
tered laser light at a fixed scattering angle can used to 
characterize the temporal behavior of the pattern am- 
plitude. The resulting trajectories have been analyzed 
quantitatively and their statistical properties have been 
extracted. 

The statistical analysis confirms that the distribution 
density of laminar phase durations r is in full agreement 
with theoretical predictions. In particular, the t~ 3 / 2 
power law describes the statistics of laminar phase dura- 
tions at the stability threshold in the conduction regime 
over four decades in r. 

The distribution density of pattern amplitudes A in the 
vicinity of the instability threshold is also in quantitative 
agreement with the predicted power law. Deviations are 
found in the experiment in the limits of low and high 
pattern amplitudes where additive noise (small A) and 
nonlinearities in the dynamic equations (large A) influ- 
ence the dynamical behaviour of the system variables. 
With increasing voltage the exponent — 1 + r\ increases. 
Theory predicts a relation rj oc X for on-off intermittence, 
which allows us to define a critical voltage U® xp from the 
slopes of the amplitude distribution functions. The crit- 
ical voltage obtained with this method agrees with the 
value derived from the analysis of the laminar phase du- 
rations, i.e. the two definitions characterize consistently 
the experimental system as well as the numerically sim- 
ulations. This provides, irrespective of the fact that the 
Lyapunov exponent is not directly accessible in the ex- 
periment, a quantitative criterion for the stability thresh- 
old U° xp of stochastically driven EHC patterns. From a 
practical aspect, the distribution density of the pattern 
amplitudes is the more sensitive measure for the deter- 



mination of U° xp . 

In the simulations of the corresponding model system, 
trajectories which are characterized by the t~ 3 / 2 and 
A^ 1 power laws are obtained when the driving parame- 
ters are chosen such that the Lyapunov exponent is zero 
(sample stability threshold |^|). This result coincides 
with theoretical predictions of the universal behaviour of 
on-off intermittency Q]. In the vicinity of the critical 
voltage the linear dependence of the exponent — 1 + r\ 
on A is also in quantitative agreement with the numeri- 
cal simulation. In the experiment, we cannot relate the 
exponent of the amplitude distribution to a Lyapunov 
exponent. However, the functional dependence of r\ on 
the reduced voltage e is in satisfactory agreement with 
the calculated data (see Fig. |l5| ). 

We note that although the statistical characterization 
of experimental and simulated data are in good quantita- 
tive agreement, the experimental and calculated thresh- 
old voltages and wave numbers can differ on an abso- 
lute scale by roughly 10% (see Fig. ||). This is mainly 
the consequence of the simplified assumptions of director 
and flow modes in the model, it is not relevant for the 
description of on-off-intermittent behaviour. 

In the power spectrum of the experimental trajecto- 
ries, a o> -1 / 2 dependence is indicated in a small frequency 
range, in agreement with predictions of general theories 
pl[ p2|, in the high frequency tail, the PSD adopts an 
ujr 2 behavior (Fig. |l6| ). In the numerical simulations 
with boundaries to the system variables, (Fig. |l7j ) we 
did not find the lo^ 1 / 2 dependence, it is reproduced only 
when boundaries of the trajectory are disregarded and 
a simulation with quasi unlimited dynamical amplitude 
range is performed. This discrepancy leads us to the con- 
clusion that the PSD is particularly sensitive to additive 
noise and the full nonlinear dynamical equations. The 
appearance of the w -1 / 2 range in the experimental data 
in apparent agreement with general predictions should 
therefore not be overestimated. 

We note that although the investigated experimental 
system represents a spatially extended dissipative sys- 
tem, it has been shown in this study that its fundamental 
statistical properties can be reproduced in simulations of 
a 2x2 evolution matrix model with global pattern am- 
plitude, thus neglecting spatial details of the pattern. In 
the case of a global driving parameter, the spatial phase 
does not play a role as long as other noise sources are 
excluded. Additive noise, however, may introduce phase 
drifts in the system |7(| [n] and is responsible for a com- 
plex spatio-temporal characteristics. A detailed spatio- 
temporal description of the system represents an ongoing 
interesting task. 

Thanks is due to Hirokazu Fujisaka for stimulating dis- 
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TABLE I: Material parameters (cgs) used in the simulations. 
The parameter set for the simulations of stochastic trajecto- 
ries is taken from the fit of threshold voltage and wave number 
characteristics for periodic AC driving, see also Fig. ^. Ex- 
perimentally data for Mischung 5 taken from |49J, BQ] 




FIG. 1: Schematic drawing of the convection rolls and director 
field in a nematic sandwich cell. A snapshot of the spatial 
modulations of director and charge fields (tp, q) in the cell 
midplane is sketched. 
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T~r — r 




Frequency zv [Hz] Frequency 1/ [Hz] 



FIG. 2: Threshold voltages and wave numbers for driving 
with periodic and stochastic (DMP) square waves of (mean) 
frequency v. Periodic excitation: experiment (o), theory 

( ), stochastic excitation: experiment (•), theory ( — ). 

The method for the experimental determination of stochastic 
thresholds is explained in the text. The symbol o indicates 
the frequency, where the stochastic measurements presented 
below have been performed. 




40 60 80 ion ]?0 140 
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FIG. 3: Transmission microscope images of the noise driven 
pattern during a burst (a) and space-time plot of a cross sec- 
tion along the wave vector (b). The horizontal axis gives 
the spatial coordinate along the pattern wave vector (direc- 
tor easy axis). Only in adjacent bursts, the spatial phase of 
the pattern appears correlated, long laminar phases destroy 
such correlations. 
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FIG. 4: Excitation waveforms and corresponding power spec- 
tral densities for periodic square wave (a) and DMP (b). 
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FIG. 5: Sketch of the experimental setup. 
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Scattering angfe 6 [rad] 

FIG. 6: Scattering image at a = 2v — 160 s _1 (jumps per sec) 
DMP excitation recorded with CCD camera, (a) 2D snap- 
shot during an intensive burst at t = 500 s, (b) scattering 
angle-time plot of intensity profiles taken at y = 0, presented 
in inverse grey scale. The constant angles of the individual 
reflexes in subsequent bursts reflects the fast and stable wave- 
length selection mechanism. An arrow marks the position of 
the photo diode, set to the most intense scattering reflex of 
the most instable wave number. 
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FIG. 7: Analytically calculated Lyapunov exponent A of the 
most instable mode, k c = 0.204 fim -1 , for DMP excitation 
with v = 80 s _1 . Material parameters are taken from the fit 
of periodic excitation thresholds and wavelengths. The value 
A = defines the critical voltage U* h = 14.2 V. 
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FIG. 8: Details of a simulated trajectory of tp and q at 
v — 80 s _1 , at A = ('conductive' regime) for a constant 
lower boundary i/) m m = 5 x 10 -3 , see Eq. (|l3|). In the low fre- 
quency regime the slow variable ip(t) keeps its sign, whereas 
q(t) oscillates synchronously with the applied field E±. 
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FIG. 9: Two simulated trajectories at A < for different 
constant background levels (a) V'mtn = 5 x 10 3 and (b) tp mill — 
5 x 10 4 . The trajectories appear to be significantly different, 
whereas the statistical analysis produce the same fundamental 
power laws. 
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FIG. 10: Same as in Fig. |8| but with background additive 
noise in tp(t). With additive noise, can change its sign 
occasionally. The pattern amplitude as well as the optical 
scattering intensities are, however, insensitive to the sign of 
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FIG. 11: Trajectories recorded with a photo diode positioned 
at the most intense second order reflex and 24 bit ADC, the 
excitation is identical with that in Fig. ^| (a) At low voltage 
(U=12.2 V) the measured raw signal A lam fluctuates around 
a background level Ao, interrupted by infrequent bursts, (b) 
Corrected intensity A(t) = \A Ta , w (t) - Ao\ at U=12.9 V where 
we assume A = (see statistical analysis). The up-down sym- 
metry of the curve is recognized. The experimental dynamic 
range is limited by saturation for large amplitudes and by 
background noise for small amplitudes, (c) At U=13.6 V the 
on-state is dominant, intermitted infrequently by breakdowns 
to the quiescent off-state. 
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Laminar Duration 2ar 

FIG. 12: Normalized distribution p(2ar) of durations r of 
laminar phases from experimental (left) and simulated (right) 
trajectories, voltages as in Figs. |ll|(a-c). Here the trajectories 
were recorded in 12 bit resolution with a 1 ms sampling rate. 
The time axis was scaled with the jump rate a = 2u. The 
threshold A on is set to 0.05 V. The t -3 / 2 power law holds 
over several decades, best agreement is found at U = U™ p 
which we assign to A = 0. 
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FIG. 13: Distribution density of the burst amplitudes A(t) 
measured over 4200 sec for same DMP voltages as in 
Figs. |ll](a)-(c). (d) for U=11.8 V (large negative A and (e) 
background at zero excitation voltage. In their middle parts, 
the distributions fit to a power laws A~ 1+v . From the curve 
with slope -1 we find the experimental voltage corresponding 
to A = . 
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FIG. 14: Amplitude distribution densities p(A t ) of the sim- 
ulated trajectories. The power law holds in the complete 
dynamical range, the fixed boundaries (Eq. (^)) generate 
abrupt edges in the support of p(A th ). The critical voltage 
Uc h = 14.2 V deviates somewhat from the experimental value. 
The graphs represent U=14.2 V (solid), 14.9 V (dotted) and 
13.5 V (dashed). 




FIG. 15: Voltage dependence of the exponent — 1 + 77 (o 
experiment, simulation). The empirical error bars illus- 
trate the uncertainty of the fit and variations between indi- 
vidual runs of the measurements. For e < 0, burst appear 
infrequently and the statistic is rather poor. A scattering ef- 
ficiency according to Eq. (^) has been used to relate measured 
optical data to the simulations (A tx \<p\ <x A 1 ^ 4 ) and assumed 
the dash-dotted line the Lyapunov exponent is zero and we 
correspond this voltage to e=0. Bigger symbols indicate the 
applied voltages depict in Fig. uA 
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FIG. 16: Power spectral density of trajectories (a)-(c) in Fig. 
|l2} AD-converter with 1 ms sample time and 12 bit resolution 
was used. General predictions for the PSD in on-off intermit- 
tency are an exponent zero for very low frequencies, (-1/2) 
for medium frequencies and (-2) for high frequencies, see Eq. 
([[til). Near the critical voltage Uc Xp , the experimental data 
indicate such a behaviour. The curve has been smoothed by 
averaging the density over intervals proportional to frequency. 
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FIG. 17: Simulated power spectrum at U* with different 
lower bounds ( Vmin = 5 x 10" 3 (•), 1CT 10 (A), 1CT 100 (o)), 
upper limit t/wx = 0.5. The theoretical exponent -1/2 is 
found only for unrealistically low background. 



